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Abstract 

Polymer  electrolyte  membrane  (PEM)  fuel  cell  performance  is  largely  controlled  by  the  microstructure  of  the  catalyst  layer.  Hence,  a  need 
exists  to  develop  fundamental  understanding  of  the  effect  of  catalyst  layer  microstructure  on  a  PEM  fuel  cell  operational  characteristics. 
Significant  obstacles  to  commercialization  of  a  PEM  fuel  cell  have  been  attributed  to  high  cost  of  the  platinum  catalyst,  slow  kinetics  of 
the  oxygen  reduction  reaction  (ORR),  and  low  catalyst  utilization.  In  order  to  address  these  limiting  factors,  a  model  has  been  developed  to 
investigate  the  effect  of  catalyst  particle  size,  loading,  utilization  and  distribution  of  catalyst  layer  components.  The  model  predicts  a  large 
increase  in  the  exchange  current  and  in  cell  potential  following  a  corresponding  reduction  in  activation  polarization  as  catalyst  particle 
size  decreases  to  a  few  nanometers.  Model  predictions  compare  well  with  published  experimental  data.  The  model  also  predicts  an  upper 
statistical  limit  of  22%  to  the  amount  of  catalyst  that  can  be  utilized  for  current  generation,  even  for  best-prepared  catalyst  layer.  This 
explains  the  observed  low  catalyst  utilization  (<25%)  reported  in  the  literature.  In  the  proposed  model,  this  statistical  limit  is  attributed  to 
the  intrinsic  characteristics  of  “randomly”  distributed  catalyst  layer  components.  To  achieve  higher  utilization,  catalyst  layer  components 
must  be  distributed  in  an  engineered  design  that  ensures  maximum  number  of  active  sites. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

High  cost  of  the  Pt  catalyst  and  the  slow  kinetics  of  the 
O2  reduction  reaction  process  necessitated  the  search  for 
alternative  ways  to  concentrate  the  Pt  catalyst  in  the  cat¬ 
alyst  layer  of  polymer  electrolyte  membrane  (PEM)  fuel 
cells.  An  effective  method  to  lower  the  Pt  loading  without 
sacrificing  performance  is  during  the  catalyst  particle  size 
to  the  nanoscale  range,  which  effectively  increases  catalyst 
surface  area  available  for  current  generation. 

In  recent  years,  nanocrystalline  materials  were  the  sub¬ 
ject  of  extensive  research  and  various  techniques  were 
developed  to  synthesize  them  [1-5].  These  materials  are 
polycrystalline  in  nature,  but  have  an  ultra-fine  grain  size 
of  the  order  of  1-100  nm.  In  such  materials  the  volume 
fraction  of  the  grain  boundaries  becomes  comparable  to 
the  volume  fraction  of  the  crystals  themselves.  As  many  as 
50%  of  the  atoms  are  located  at  interfacial  boundaries.  Due 
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to  this  large  volume  fraction  of  grain  boundaries,  and  the 
high  surface  area  of  these  materials,  they  exhibit  different 
properties  than  conventional  engineering  alloys  [6]. 

Interest  in  nano  structured  materials  as  potential  catalysts 
is  evident.  The  Literature  shows  that  nano  structured  cata¬ 
lysts  exhibit  improved  catalytic  activities  [7-13].  However, 
because  the  processing  and  the  application  of  nanostructures 
to  catalysts  is  in  its  infancy,  more  work  is  needed  to  provide 
a  fundamental  understanding  of  this  new  class  of  catalysts. 
Along  with  experimental  observation,  theoretical  modeling 
is  needed  to  fully  exploit  microstructural  parameters  in  an 
attempt  to  optimize  fuel  cell  performance. 

Several  mathematical  models  have  been  proposed  in  the 
literature  for  PEM  fuel  cells;  these  models  have  concentrated 
on  various  aspects  of  PEM  fuel  cell  operation,  and  their 
effect  on  cell  performance.  Models  developed  considered 
effects  such  as  mass  transport,  porosity,  gas  pressures,  water 
management,  reaction  kinetics,  etc.  A  brief  review  of  several 
of  these  models  is  given  below. 

Mann  et  al.  [14]  and  Amphlett  et  al.  [15]  developed 
a  parametric  model  of  a  single  PEM  cell  by  using  a 
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mechanistic  approach,  and  a  number  of  grouped  parameters 
were  identified  and  fit  to  an  empirical  data  measured  from  a 
Ballard  Mark  IV  single  cell.  Similarly,  Kim  et  al.  [16]  used 
an  empirical  equation  to  describe  the  performance  data  of  a 
PEM  fuel  cell  over  a  complete  operating  range.  Grot  et  al. 
[17]  and  Fuller  and  Newman  [18]  incorporated  strategies  for 
thermal  and  water  management.  Marr  and  Li  [19]  extended 
Grot  et  al.  model,  and  investigated  catalyst  utilization  as  well 
as  the  optimal  composition  and  structure  for  the  cathode  cat¬ 
alyst  layer;  such  as,  catalyst  loading,  catalyst  type,  catalyst 
layer  thickness,  void  fraction,  and  ionomer  content.  Springer 
et  al.  [20]  developed  an  isothermal,  one-dimensional,  steady 
state  model  for  water  transport  through  a  complete  PEM 
fuel  cell  based  on  experimental  results;  their  modeling 
results  provided  useful  insight  into  the  cell’s  water  trans¬ 
port  mechanisms  and  their  effect  on  the  cell  performance. 
Bernardi  and  Verbrugge  [21,22]  developed  a  comprehensive 
mathematical  model  for  a  PEM  fuel  cell  from  fundamental 
transport  properties.  Nguyen  and  White  [23]  formulated  a 
quasi-two-dimensional  model  to  account  for  the  variations 
of  mass  and  heat  transfer  between  the  electrode  and  reactant 
gas  mixture  in  the  flow  channel.  Recently,  Weisbrod  et  al. 
[24]  presented  a  simplified  through-the-electrode  model, 
including  water  saturation  and  transport  within  the  elec¬ 
trode  and  reaction  kinetics  within  the  cathode  catalyst  layer. 
Work  of  Lee  et  al.  [25]  presents  the  most  comprehensive 
form  of  an  empirical  model  produced  to  date  to  predict 
the  current,  voltage  relationship  for  the  typical  PEM  fuel 
cell.  Gottesfeld  and  Zawodzinski  [26]  have  also  achieved 
considerable  success  in  PEM  fuel  cell  modeling  over  the 
last  decade,  resulting  in  increasingly  complex  predictors 
of  cell  performance  based  on  interdependencies  of  process 
variables  and  the  nature  of  transport  processes.  Eikerling 
and  Kornyshev  [27]  proposed  a  phenomenological  model 
that  investigated  the  effect  of  O2  diffusion,  proton  con¬ 
duction  and  reaction  kinetics  in  the  catalyst  layer.  You 
and  Liu  [28]  presented  a  pseudo-homogeneous  model  for 
the  cathode  catalyst  performance  in  PEM  fuel  cell  derived 
from  the  basic  mass  current  balance  by  the  control  volume 
approach. 

To  date,  a  comprehensive  model  that  directly  relates 
catalyst  layer  microstructural  refinement  and  catalyst  uti¬ 
lization  to  PEM  fuel  cell  performance  is  non-existent.  The 
objective  of  this  work  is  to  develop  a  mathematical  model 
relating  catalyst  structure  refinement,  loading,  utilization 
and  distribution  to  polarization  effects  and  overall  PEM 
fuel  cell  performance.  In  this  paper,  a  modified  expression 
for  the  exchange  current  in  the  Butler- Volmer  equation  is 
developed.  The  new  expression  incorporates  a  geometrical 
factor  which  includes  both  catalyst  particle  size  and  load¬ 
ing  for  a  random  distribution  of  catalyst  layer  components. 
Furthermore,  a  new  parameter  “utilization  coefficient”  has 
been  introduced,  which  represents  the  fraction  of  catalyst 
surface  area  located  in  electrocatalytically  active  sites.  The 
discussion  in  this  paper  may  help  in  the  optimization  and 
design  of  cathode  catalyst  layer. 


2.  Model  description 

In  this  study  the  catalyst  layer  is  assumed  to  consist  of  a 
mixture  of  catalyst  platinum,  ionomer  membrane  electrolyte, 
carbon  support  and  void  space  (network  of  pores).  Although 
the  catalyst  layer  is  relatively  small,  it  is  the  heart  of  the 
fuel  cell.  Here  fuel  and  oxidant  react  electrochemically  to 
produce  electrical  energy.  Four  different  media  are  present 
for  the  function  of  a  catalyst  layer:  an  interconnected  solid 
phase  of  the  carbon,  which  serves  as  the  supplying  network 
for  the  electrons;  a  network  of  gas  pores,  which  supplies 
oxygen  to  the  reaction  sites;  an  electrolyte  network  provides 
the  pathways  for  the  protons  to  the  catalyst  places;  and  a 
catalytic  surface  for  electrochemical  reaction  to  take  place 
(Fig.  1). 

The  performance  of  H2/O2  (fuel/oxidant)  PEM  fuel  cell 
is  limited  primarily  by  the  slow  rate  of  the  O2  reduction  half 
reaction,  O2  +  4H+  +  4e_  — >►  2H2O,  which  is  many  more 
times  slower  than  the  H2  oxidation,  2H2  — >  4H+  +  4e_ . 
When  current  flows  through  a  fuel  cell,  the  potential  drop 
due  to  polarization  effect  which  is  primarily  controlled  by 
the  oxygen  reduction  reaction  (ORR)  rate. 

In  PEM  fuel  cell,  the  three-phase  interface  between  the 
electrolyte,  the  electrode  and  the  reactant  gas  (within  a  per¬ 
colating  network  of  ionomer,  oxygen  and  electronic  con¬ 
ducting  medium)  is  important.  Only  catalyst  located  in  this 
region  is  electrochemically  active,  i.e.  capable  of  generating 
current.  Clearly,  the  rate  at  which  this  happens  will  be  pro¬ 
portional  to  the  area  of  the  electrode.  Therefore,  increasing 
the  catalyst  surface  area  by  reducing  catalyst  particle  size 
will  give  rise  to  increased  reaction  rates. 

Previous  models  assume  that  the  half-cell  reaction  has 
equal  probability  to  occur  anywhere  on  the  electrode  sur¬ 
face  (apparent  area).  In  PEM  fuel  cell,  real  area  apparent 
area  and  not  all  electrode  real  area  is  available  for  reaction 
but  only  active  sites  are  utilized  for  current  generation.  The 
exchange  current,  which  reflects  the  activities  of  electrode, 
depends  largely  on  the  active  surface  area  of  the  electrode. 
In  this  model,  the  exchange  current  in  the  Butler- Volmer 
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Fig.  1.  Schematic  diagram  of  an  MEA  illustrating  the  catalyst  layer 
microstructure. 
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equation  has  been  modified  to  include  active  surface  area, 
which  accounts  for  catalyst  particle  size,  loading  and  utiliza¬ 
tion  to  give  improved  description  of  the  electrode  reaction 
kinetics. 


3.  Model  formulation 


the  entire  real  or  nominal  electrode  surface  areas.  The  net 
reaction  rate  is,  therefore,  the  difference  between  forward 
and  backward  reaction  rates,  i.e., 

-» =•"-»-= k<c*s  -  ‘-c-s = ( ilk)  ~  (ss)  <5> 

rearranging, 


Fuel  cell  potential  as  a  function  of  polarization  effects 
[29]  for  a  single  cell  is  given  by: 

4x11  =  ^eq  T  ?7c,act  T-  ^?a,act  T  Y]c  ,conc  T  0  a,  cone  T"  Johnnie  (1) 

where  £Ceii  is  the  irreversible  measured  voltage  in  the  cell 
(when  a  net  current  passes  through  the  circuit),  EQ q  is  the 
Nernst  equilibrium  potential  for  the  cell  or  reversible  ther¬ 
modynamic  potential  also  known  as  open  circuit  voltage 
(OCV)  (when  no  net  current  passes  through  circuit).  r]c, act 
and  77a, act  are  cathode  and  anode  activation  overvoltages 
associated  with  kinetics  limitations  to  the  charge  transfer 
process  at  the  electrodes;  0c,conc  and  o  a, cone  are  concentra¬ 
tion  overvoltages  resulting  from  depletion  of  reactants  in 
the  vicinity  of  the  electrodes  due  to  slow  diffusion  from 
bulk  solution;  and  770hmic  is  the  resistive  overvoltage  asso¬ 
ciated  with  the  internal  resistance  in  electrode,  electrolyte 
and  collectors.  Experimental  data  to  date  on  PEM  fuel 
cells  shows  that  the  anode  activation  polarization  is  neg¬ 
ligible  with  respect  to  the  cathode  activation  polarization 
due  to  slow  kinetics  of  the  oxygen  reduction  reaction  at  the 
cathode. 

Consider  the  following  hypothetical  half-cell  electrode  re¬ 
action  (cathode  reaction): 

R  ne  P  (2) 


The  forward  and  backward  reactions  rates  Vf  and  Vb 
(mol/cm2  s),  respectively,  are  as  follows: 


Vf  =  kfCRS  = 


nFA 


Vb  =  kbCpS  = 


Jb_ 

nFA 


where  kf  and  kb  are  the  rate  constants  for  the  forward  and  the 
backward  reactions  in  (cm-1  s-1);  Cr  and  Cp  are  reactant 
and  product  concentrations  (mol/cm3),  respectively.  In  this 
model  a  parameter,  S,  is  incorporated  into  the  reaction  rate 
expression,  a  “geometrical  factor”,  that  is  a  function  of  the 
catalyst  layer  microstructure  (catalyst  particle  size,  loading 
and  utilization)  and  is  equal  to  the  total  “active”  surface  area 
(cm2).  If  and  4  are  the  forward  (cathodic)  and  backward 
(anodic)  current  contributions  to  the  total  half-cell  current 
in  Amperes  or  Coulombs/s.  A  is  the  electrode  nominal  sur¬ 
face  area  (cm2),  F  is  Faraday  constant  (96,485  C/mol)  and 
n  is  the  stoichiometric  number  of  electrons  consumed  in 
the  electrode  reaction.  The  need  for  a  geometrical  factor 
comes  from  the  fact  that,  for  PEM  fuel  cell,  reaction  takes 
place  only  at  active  sites  within  the  catalyst  layer  and  not  at 


4et  =  If  —  h  =  nFAS[kfCR  —  hCp] 


where  4e t  is  the  net  half-cell  current.  Writing  the  rate  con¬ 
stants  as  function  of  the  standard  rate  constant  ( k° )  and  the 
half-cell  overpotential  (E  —  Eeq), 


kf  =  k°  exp 


—anF(E  -  Eeq) 
RT 


kb  ~  k°  exp 


a)nF(E 

rt~ 


where  of  is  a  charge  transfer  coefficient  and  E  is  the  half-cell 
potential.  Substituting  Eqs.  (7)  and  (8)  into  Eq.  (6)  and 
assuming  no  mass  transfer  effects  (good  convection),  this 
yields  a  “modified”  Butler- Volmer  equation: 


I  =  L 


o 


(  —anFr]\ 

pyn  1  I  pvn  „ 

e 

1 

H 

exP  1  J  exP 

RT 

where  77  is  the  activation  overpotential  and  is  equal  to  E  — 
£eq,  and  I0  is  a  modified  value  for  the  exchange  current: 

/0  =  nFA  Sic  C*' 1 ~a)  Cjf  (10) 


where  and  Cp  the  are  bulk  concentrations  of  reactants 
and  products,  respectively.  For  the  cathodic  reaction  where  77 
is  relatively  large,  the  second  term  in  Eq.  (9)  becomes  negli¬ 
gible.  Hence,  the  modified  Butler- Volmer  equation  reduces 
to: 


(11) 


Dividing  by  the  nominal  area  A  and  rearranging,  then  the 
activation  polarization  for  the  oxygen  reduction  reaction  at 
the  cathode  will  become: 


(12) 


where 

io.c  =  nFSk°C*f-ac)C*f-ac)C^0  (13) 

Eq.  (12)  is  known  as  Tafel  equation  and  the  term 
( RT/nacF )  is  the  Tafel  slope.  ac  represents  the  cathodic 
transfer  coefficient,  Cq2,  C^+  and  C^oQ  are  the  bulk  con¬ 
centrations  of  O2,  H+  and  H2O  in  the  PEM  fuel  cell, 
respectively,  i  is  current  density  and  z0,c  is  the  modified 
exchange  current  density  for  the  cathode  reaction  (oxygen 
reduction  reaction  on  Pt  catalyst). 
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3.1.  Determination  of  the  geometrical  factor 


Assuming  that  the  catalyst  layer  consists  of  nc  randomly 
distributed  spherical  catalyst  particle  and  that  the  single  cat¬ 
alyst  particle  surface  area  is  given  by  Sc  (cm2),  then  the  total 
active  surface  area,  S,  is  given  by: 

S  =  yncSc  (14) 


where  y  is  a  “utilization  coefficient”  which  is  the  fraction 
of  catalyst  surface  satisfying  active  area  requirements  for 
current  generation  (three-phase  contact).  That  is,  particles 
which  are  not  in  contact  with  the  electronic  conducting  phase 
(catalyst  and  carbon  support),  electrolyte  and/or  gas  perco¬ 
lating  networks  will  not  contribute  to  the  active  surface  area 
and  no  current  is  generated  from  these  areas.  Substituting 
the  surface  area  of  a  spherical  particle  of  radius,  rc  (cm), 
into  Eq.  (14),  we  get 

S  =  4n  r^ncy  (15) 

and  the  total  volume  of  catalyst  particles,  Vc  (cm3),  is  given 
by: 

mc  /4\  o 

i,c= * = (5)  ^  (i6> 


where  mc  (g)  and  pc  (g/cm3)  are  the  catalyst  mass  and  den¬ 
sity,  respectively.  Substituting  nc  from  Eq.  (16)  into  Eq.  (15), 
the  active  surface  area  becomes: 


3  ymc 
rcpc 


(17) 


Now,  substituting  Eqs.  (13)  and  (17)  into  Eq.  (12): 


he, act 


In 


3ymcnFk° Cq)  “c)C*(1  ac) 


H+ 


/-<*ac 

ch2o 


ircpc 


(18) 


3.2.  Determination  of  the  “Utilization  Coefficient” 

Consider  three  sets  of  balls  each  has  a  different  color, 
where  nw ,  n\>  and  nY  are  the  numbers  of  white,  black  and 
red  balls,  respectively.  If  you  were  to  pick  three  balls  from  a 
randomly  mixed  set  of  nw  +  n\>  +  nv  balls,  then  the  question 
is:  what  is  the  probability  of  picking  one  ball  from  each 
color?  In  contrast  to  our  catalyst  layer  structure,  where  colors 
represent  electrolyte,  catalyst  and  gas  regions,  the  question 
would  become:  what  is  the  probability  of  finding  the  three 
regions  in  contact? 

Therefore,  a  utilization  coefficient  (y),  in  statistical  terms, 
can  be  defined  as  the  probability  of  having  contact  between 
catalyst  (supported  on  an  electronic  conducting  carbon  net¬ 
work),  electrolyte  and  gas  in  a  randomly  distributed  con¬ 
stituents  within  the  catalyst  layer.  Then  y  can  be  calculated 
as  follows: 


n  1 


’n2  Cr2 
U[  Crt 


’»3 


(19) 


where  nCr  is  the  number  of  ways  in  which  r  objects  can 
be  selected  from  a  set  of  n  distinct  objects,  also  called  the 
number  of  “combinations”  and  „tC n  is  the  total  number  of 
combinations.  For  a  catalyst  layer  Eq.  (19)  can  be  written  as: 


U[ 


(20) 


where  UcC  1  represents  the  number  of  combinations  of  find¬ 
ing  one  catalyst  particle  from  a  set  of  nc  catalyst  particles, 
neC\  and  nsC\  are  defined  in  a  similar  manner.  „tC3,  gives 
the  total  number  of  combinations  of  finding  three  regions 
from  a  set  of  nt  =  n&  +  nc  +  ns.  Although  there  are  struc¬ 
tural  differences  between  the  catalyst  and  the  ionomer,  it  is 
reasonable  to  assume  that  at  the  contact  zone  electrolyte  can 
be  regarded  as  “particle-like”  medium,  y  given  by  Eq.  (20) 
varies  between  0  and  1  and  represents  the  fraction  of  cata¬ 
lyst  surface  area  that  lies  in  the  three-phase  region.  It  should 
be  emphasized  that  the  above  argument  is  valid  only  in  a  re¬ 
gion  above  the  percolation  threshold.  That  is,  assuming  that 
the  active  sites  (three-phase  regions)  have  access  to  the  elec¬ 
tronic  conduction  path,  oxygen  pathway  and  an  electrolyte 
network.  Replacing  y  in  Eq.  (18)  by  its  value  from  Eq.  (20), 
the  modified  activation  polarization  becomes: 


0  c, act 


In 


U[  c3 


X 


m  nFk°C*^  occ)  ^ac 

mcnrn  tQ2  <-H+  <-h2Q 

Ire  Pc 


(21) 


Using  analysis  in  reference  [15],  it  can  be  shown  that  the 
concentration  overpotential  is  given  by: 


7?  cone  — 


(22) 


where  i\  known  as  the  limiting  current  density  for  the  elec¬ 
trode  that  has  the  lowest  limiting  current  density.  The  lim¬ 
iting  current  density  is  the  current  at  which  the  fuel  is  used 
up  at  a  rate  equal  to  its  maximum  supply  speed.  For  the 
Ohmic  polarization,  potential  drop  is  simply  proportional  to 
the  current  and  is  given  by: 


0  ohmic  —  lr  (23) 

where  r  is  the  area  specific  resistance  in  the  cell  (£2  cm2). 

Substituting  all  polarization  effects  (Eqs.  (21)-(23))  into 
Eq.  (1), 
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Ece\\  —  Eq  q  —  (/  +  in)r  + 


RT 


x  In 


nacE 
[npC  1  ^1  ^l] 


X 


RT\ 

+  I  —  I  In 


nt^3 
-ac 

02  ^H+ 


m  nFk°C ac)^*a;c 
rncnrt c  e  n.  ^u+  ^H20 


(/  +  G)rcpc 

0  +  in) 


nF 


1  - 


(24) 


where  /„  is  a  term  added  to  account  for  the  internal  and  fuel 
cross  over  equivalent-current  density  [30].  Eq.  (24)  gives 
the  measured  cell  voltage  as  a  function  of  catalyst  particle 
size,  loading  and  utilization  for  a  randomly  distributed  con¬ 
stituents  in  the  catalyst  layer.  Feq  can  be  calculated  from  the 
molar  free  energy  of  formation  of  H2O,  AGf,  using  the  rela¬ 
tionship  Eq  q  =  —  AGf/ nF,  where  n  is  the  number  of  electrons 
per  each  molecule  of  fuel  consumed,  this  can  also  be  deter¬ 
mined  from  Nernst  equation  as  a  function  of  concentrations. 


4.  Results  and  discussion 


The  influence  of  different  micro  structural  parameters  in  a 
catalyst  layer  can  be  studied,  in  depth,  by  examining  model 
predictions  given  by  Eq.  (24).  Substituting  typical  values 
[15,27,30-32]  of  the  parameters  in  Eq.  (24)  for  Pt  catalyst: 


£<*11  =  1.2-  (i  + 2)3  x  1(T5 


+  0.06  In 
+  0.05  In 


9.37  x  10  3ymc 


rc(i  +  2) 
"/  +  2" 
900 


(25) 


where  units  for  ECq\\,  i ,  mc  and  rc  are  volts,  mA/cm2,  g  and 
cm,  respectively;  y  is  given  by  Eq.  (20).  A  peak  value  of  y 
achieved,  for  a  random  mix,  when  the  number  of  catalyst, 
gas  and  electrolyte  regions  are  equal.  Fig.  2  gives  the  value 
of  y  as  a  function  of  the  number  of  gas  and  electrolyte 
regions,  for  a  constant  value  nc  of  103.  It  can  be  seen  from 
this  figure  that  when  nQ  =  ns  =  nc  a  spike  appears  which 
represents  the  highest  y  value.  Furthermore,  for  nQ  =  ns 
=  nc  =  n.  Fig.  3  reveals  that  for  a  large  n  (>100),  the  highest 
value  of  y  achievable  is  approximately  0.22.  In  other  words, 
according  to  this  model,  only  22%  or  less  of  catalyst  surface 
is  available  for  current  generation,  since  only  particles  in 
the  three-phase  region  are  electrochemically  active.  In  the 
following  discussion,  a  value  of  y  =  0.22,  which  is  the  most 
optimistic  value  of  the  “utilization  coefficient”,  is  used. 

A  crucial  factor  in  improving  fuel  cell  performance  is 
therefore  to  increase  the  value  of  the  exchange  current,  espe¬ 
cially  at  the  cathode.  According  to  the  model,  the  exchange 
current  is  highly  influenced  by  catalyst  particle  size  and 
loading.  The  log/log  plot  in  Fig.  4  demonstrates  the  effect  of 
catalyst  particle  size  and  loading  on  the  exchange  current. 
As  the  particle  size  is  reduced  from  1  p,m  to  1  nm,  i0  in¬ 
creases  by  three  orders  of  magnitude.  The  effect  is  due  to  a 
rapid  increase  in  catalyst  surface  area,  in  particular,  when  rc 
drops  below  40  nm  as  illustrated  in  Fig.  5.  Furthermore,  for 
a  given  i0,  loading  can  be  lowered  from  100  to  0.1  mg/cm2 
provided  that  particle  size  is  reduced  from  0.5  mm  to  15  nm. 
That  is,  lowering  catalyst  loading  can  be  compensated  for 
by  using  smaller  catalyst  particle  size  without  effecting  fuel 
cell  performance. 

The  i-v  characteristics  in  Fig.  6  reveal  the  dependence  of 
fuel  cell  potential  on  catalyst  particle  size.  The  plot  clearly 
shows  how  reducing  catalyst  particle  size  reduces  overpo¬ 
tential  and  increases  fuel  cell  voltage.  Cell  voltage  rises 
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Fig.  2.  Catalyst  utilization  coefficient  as  a  function  of  ne  and  ns  for  nc  =  103. 
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n 

Fig.  3.  Catalyst  utilization  coefficient  vs.  n,  where  n  =  nc  =  ne  =  ns. 


from  0.16  to  1.00  V  as  rc  decreases  from  100  fxm  to  1  nm 
at  300mA/cm2.  A  similar  behavior  is  observed  when  ex¬ 
amining  the  effect  of  loading  on  cell  performance  (Fig.  7). 
As  expected,  for  a  given  particle  size,  increasing  catalyst 
loading  has  the  effect  of  shifting  the  i-v  plot  towards  higher 
potential.  However,  according  to  the  model,  higher  cell 
specific  power  can  be  achieved  when  using  lower  catalyst 
loading  (Fig.  8).  The  effect  of  catalyst  particle  size  on  % 
power  loss  (Fig.  9)  shows  significant  increase  in  %  poten¬ 
tial  loss  accompanying  the  increase  in  catalyst  particle  size. 
Typical  PEM  fuel  cell  has  loading  of  0.5  mg/cm2  and  about 
45%  loss  in  potential  at  a  current  density  of  300mA/cm2. 
It  can  be  seen  that,  for  the  same  potential  loss  of  45%, 
catalyst  loading  of  0.01  mg/cm2  is  possible  if  particle  size 
is  reduced  to  2  nm,  which  is  a  large  reduction  in  Pt  cost. 

The  present  model  is  tested  against  reported  data  pub¬ 
lished  in  the  open  literature.  Table  1  compares  fuel  cell 


output  potentials  predicted  from  the  model  to  published 
experimental  values  [33].  It  clearly  shows  that  the  proposed 
model  is  in  good  agreement  with  experimental  data.  For 
various  catalyst  particle  sizes  and  loadings,  the  differences 
between  model  predicted  potentials  and  experimental  re¬ 
sults  are  in  the  range  of  7.1-15.9%.  The  effect  of  catalyst 
particle  size  below  1  nm  on  the  intrinsic  electrocatalytic 
activities  is  controversial  [37]  as  a  consequence  the  present 
model  is  only  valid  above  a  1  nm  particle  size. 

The  importance  of  catalyst  utilization,  represented  in  this 
model  by  the  “utilization  coefficient”,  y,  cannot  be  overem¬ 
phasized.  As  mentioned  above,  in  a  random  distribution  of 
catalyst  layer  components,  a  maximum  of  22%  of  compo¬ 
nents  (catalyst,  gas,  and  electrolyte)  have  direct  contact,  i.e., 
at  best,  only  22%  of  catalyst  is  utilized.  In  other  words,  the 
probability  of  having  a  Pt  particle  in  active  site  is  <22%. 
This  upper  “statistical”  limit  to  the  amount  of  catalyst 


Fig.  4.  The  effect  of  the  exchange  current  on  catalyst  particle  size  and  loading. 
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Fig.  5.  Particle  size  effect  on  active  surface  area  for  a  catalyst  utilization  coefficient  of  0.2  and  catalyst  loading  of  1  mg/cm2 


Catalyst  Particle  Size  |nm) 

Fig.  6.  Effect  of  Pt  catalyst  particle  size  on  the  cell  polarization  for  catalyst  loading  of  1  mg/cm2. 


Table  1 

Comparison  of  model  predictions  to  experimental  data 


Current  density  (mA/cm2) 

Cell  potential,  experiment  (volts)a 

Cell  potential,  present  model  (volts) 

Difference  (%) 

Pt  loading:  2.64 mg/cm2;  particle  size: 

3.9  nm 

600 

0.58 

0.69 

15.9 

800 

0.52 

0.56 

7.1 

Pt  loading:  3.96 mg/cm2;  particle  size: 

8.8  nm 

400 

0.74 

0.85 

12.9 

600 

0.68 

0.76 

10.5 

800 

0.58 

0.65 

10.7 

Pt  loading:  5. 12 mg/cm2;  particle  size: 

2.5  nm 

400 

0.71 

0.89 

9.2 

800 

0.61 

0.69 

11.6 

a  Reference  [33]. 
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Fig.  7.  Catalyst  loading  effect  on  the  i-v  curve  in  PEM  fuel  cell  for  catalyst  particle  size  of  1  nm. 


(1  nm)  Catalyst  Particle  Size  (cm) 


Fig.  8.  PEM  fuel  cell  specific  power  as  a  function  of  catalyst  particle  size  and  loading  at  current  density  of  300mA/cm2. 


utilized  in  a  PEM  fuel  cell  is  “intrinsic”  to  the  random  mix 
of  components  in  the  catalyst  layer,  which  explains  the  low 
catalyst  utilization  reported  in  the  literature  [34,35].  Pt  uti¬ 
lization,  even  in  the  best  electrodes,  is  quit  low  (10-25%) 
[36].  The  problem  of  low  Pt  utilization  has  been  attributed 
to  a  number  of  reasons  [27,36]:  too  high  or  too  low  Nation 
content  results  in  loss  of  connectivity  of  catalyst  parti¬ 
cles  to  membrane  and  diffusion  layer;  high  PTFE  content 
in  the  catalyst  layer  increases  hydrophobicity,  resulting 
in  improved  oxygen  supply  and  product  water  removal 
at  the  expense  of  the  Nation  percolation  which  effects 
the  supply  of  protons;  Pt  particles  located  in  very  fine 
pores  are  not  accessible  to  protons  due  to  high  molecular 
weight  of  Nation  and  that  is  not  able  to  penetrate  small 
pores. 


“On  the  other  hand,  high  catalyst  utilization  has  also  been 
reported  [38-42].  The  higher  catalyst  utilization  was  at¬ 
tributed  to  more  efficient  catalyst  layer  structures  and 
MEA  preparation  methods.  For  example,  Xiaoliang  et  al. 
[38]  found  that  the  utilization  of  catalyst  particles  in  the 
immersed  and  brushed  E-TEK  electrode  is  77.8  and  22%, 
respectively.  This  shows  that  when  the  E-TEK  is  impreg¬ 
nated  by  immersion  so  that  sufficient  proton  passageways 
are  provided  to  the  catalyst  particles,  the  platinum  utiliza¬ 
tion  can  increase  to  77.8%,  which  is  higher  than  that  of 
thin-film  catalyst  layer.” 

The  present  model  predicts  that,  even  if  optimum  con¬ 
nectivity  and  percolation  are  achieved,  the  Pt  utilization  will 
remain  low  due  to  intrinsic  characteristics  associated  with  a 
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(1  nm) 


Particle  Size  (cm) 


Fig.  9.  Effect  of  particle  size  and  loading  on  loss  in  cell  voltage  at  current  density  of  300mA/cm2. 


randomly  distributed  catalyst  layer  components.  The  utiliza¬ 
tion  coefficient  can  only  be  raised  through  designing  cata¬ 
lyst  layer  in  a  way  that  provides  higher  three-phase  contact 
than  achieved  from  “random”  conventional  mixing  of  cata¬ 
lyst  components.  Innovative  designs  and  methods  of  applica¬ 
tions  to  make  better  use  of  catalyst  are  needed.  PEM  fuel  cell 
requires  a  tailored  design  of  electrodes  for  optimal  catalyst 
placement,  hence,  it  is  expected  that  these  better-prepared 
catalyst  layer  will  have  higher  values  of  y  for  same  catalyst 
loading. 

The  effect  of  catalyst  utilization  coefficient  is  illustrated 
in  Fig.  10.  A  hypothetical  45%  increase  in  cell  output 


voltage  from  0.58  to  0.84  V  (at  i  =  300mA/cm2)  as  y  is 
raised  from  0.02  to  1.  The  model  also  allows  for  the  de¬ 
termination  of  the  utilization  coefficient  experimentally. 
A  utilization  map,  which  is  active  surface  area,  S,  ver¬ 
sus  ( mjrc )  ratio  for  different  y  values  (from  Eq.  (17)), 
is  illustrated  in  Fig.  11.  This  diagram  provides  a  valuable 
tool  that  allows  the  determination  of  y  experimentally 
using  independently  measured  S ,  mc  and  rc.  The  signif¬ 
icance  of  this  map  is  that  it  provides  a  way  to  assist  in 
the  optimization  of  catalyst  layer  by  trying  different  de¬ 
signs  and  testing  their  utilization  coefficient  using  this 
map. 


Fig.  10.  The  dependence  of  the  PEM  fuel  cell  i-v  characteristics  on  the  catalyst  utilization  coefficient. 
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Fig.  11.  Catalyst  utilization  map  illustrating  the  dependence  of  the  active  surface  area  on  the  catalyst  loading  to  particle  size  ration  for  various  catalyst 
utilization  coefficients. 


5.  Conclusions 

•  A  model  has  been  developed  to  relate  PEM  fuel  cell  per¬ 
formance  to  catalyst  particle  size,  loading  and  utilization 
based  on  randomly  distributed  catalyst  layer  components. 

•  A  modified  Butler- Volmer  equation  has  been  developed 
and  proved  to  be  in  good  agreement  with  experimental 
data. 

•  The  model  predicts  a  0.84  V  increase  (at  300mA/cm2)  in 
cell  potential  as  Pt  particle  size  is  reduces  from  100  pim 
down  to  1  nm. 

•  The  model  also  predicts  that  the  exchange  current  at  1  nm 
Pt  particle  size  is  three  orders  of  magnitude  higher  than 
that  at  1  p,m. 

•  A  utilization  coefficient  has  been  incorporated  into  the 
model,  whose  value  is  a  function  of  the  distribution  of 
components  in  the  catalyst  layer. 

•  The  model  predicts  an  upper  limit  of  22%  of  the  amount 
of  catalyst  that  can  be  utilized  in  a  randomly  distributed 
catalyst  layer  constituents. 

•  Catalyst  utilization  can  only  be  improved  through  inno¬ 
vative  designed  placement  of  catalyst  layer  components. 

•  A  utilization  map  has  been  developed  and  constructed, 
which  can  be  used  to  predict  and  design  a  catalyst  mi¬ 
crostructure  that  optimizes  fuel  cell  performance. 
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